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Conclusive mathematical arguments are presented supporting the ratchet conjecture [R. Chacon, 
J. Phys. A 40, F413 (2007)], i.e., the existence of a universal force waveform which optimally 
enhances directed transport by symmetry breaking. Specifically, such a particular waveform is shown 
to be unique for both temporal and spatial biharmonic forces, and general (non-perturbative) laws 
providing the dependence of the strength of directed transport on the force parameters are deduced 
for these forces. The theory explains previous results for a great diversity of systems subjected to 
such biharmonic forces and provides a universal quantitative criterion to optimize any application 
of the ratchet effect induced by symmetry breaking of temporal and spatial biharmonic forces. 

PACS numbers: 05.60.-k 

Symmetry principles play a fundamental role in the laws of nature by constituting a synthesis of those regularities 
that are independent of the specific dynamics. Since the beginning of scientific thinking, most of these symmetry 
principles have been deeply associated with certain "principles of economy" of nature, such as the principle of least 
action and other variational principles [1]. Curie's principle [2] is just one example of such a symmetry principle. It has 
been applied for instance to the ubiquitous phenomenon of directed transport by symmetry breaking without any net 
external force— the so-called ratchet effect [3,4,5]. Directed ratchet transport (DRT) is nowadays understood as a result 
of the interplay of nonlinearity, symmetry breaking, and non-equilibrium fluctuations, where these fluctuations may 
include temporal noise [6] , spatial disorder [7] , and quenched temporal disorder [8] . The space-time symmetries of some 
generic equations of motion, which have to be broken to allow the appearance of DRT, have recently been proposed 
[9,10,11] in accordance with Curie's principle. While this symmetry analysis provides useful necessary conditions on 
the ac (temporal) forces and the static (spatial) potential for DRT to appear, no information at all concerning its 
strength can be obtained from such a symmetry analysis. This is not altogether surprising since Curie's principle per 
se, in contrast to variational principles, lacks a mathematical formulation connecting quantitatively a measure of the 
"degree of symmetry" of causes with that of effects, thus impeding the observation of its possible connection (if any) 
with some natural "principle of economy" [12]. In other words, Curie's principle must be considered as belonging to 
the philosophical rather than the scientific realm [13]. To overcome this difficulty in the context of DRT, there has 
recently been proposed a quantitative measure of the degree of symmetry breaking (DSB) on which the strength of 
DRT must depend [14]. It should be stressed that this quantitative relationship between cause (symmetry breaking) 
and effect (DRT)— the so-called DSB mechanism [14]— is absolutely absent in the formulation of Curie's principle. 
Once one assumes that the breakage of the relevant (for each particular equation of motion) space-time symmetries 
can be quantified, and that the strength of the resulting DRT (hereafter referred to as (V)) is proportional to this 
degree of breakage, the following questions naturally arise: Can one find the parameter values of the periodic zero- 
mean forces involved that maximally break the relevant symmetries and thus optimally enhance DRT? Are these 
values universal, implying therefore the existence of a unique optimal force waveform? 

In this Letter, a positive response to these questions is provided by completing the mathematical proof of the ratchet 
conjecture [14], i.e., the existence of a universal force waveform which optimally enhances DRT. The biharmonic 
temporal force F(t')/e = i] cos (ujf + ipi) + (1 — -q) cos (2wt' + <^ 2 ) and the biharmonic spatial potential U(x)/e = 
77 sin {kx' + ipi) + (1 — T]) sin (2kx' + tp 2 ), (e > 0, r\ £ [0, 1]), have been (and still are) overwhelmingly used as standard 
models in analytical research on ratchets [15]. A simple re-scaling of these biharmonic functions yields 

f(t) = r, cos (t) + (I-77) cos (2t + ¥> e //), (1) 
g{x) = 77 cos (x) + 2 (1 — 77) cos (2x + <p e ff) , (2) 

where t = ujt' + (pi,x = kx' + <pi,f(t) = F(t)/e, g(x) = — (efc) 1 dll(x)/dx, and <p e ff — 9?2 — 2</?i is the effective 
phase. Thus, the breakage of the three relevant symmetries associated with the forces (1) and (2) (i.e., the shift 
symmetry and the two reversal symmetries, see Refs. [6,9,10,11]) is controlled by only two parameters: 77 (relative 
amplitude of the two harmonics) and feff- But, unfortunately, changing these parameters implies also changing the 
amplitude and symmetry of the positive and negative parts of the biharmonic forces with respect to the symmetric 
cases 77 = {0, 1} (see Figs. 1(a), 2(a), and 3(a)). Since the strength of any transport (induced by symmetry breaking 
or not, i.e., by non-zero- mean forces) depends upon the amplitude of the driving forces, one concludes that these two 
effects of changing 77 or <p e ff overlap, so that one will find it difficult to distinguish the contribution to transport 
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that is purely due to symmetry breaking. A first step towards clarifying this problematic situation and answering the 
above questions has been to study the breakage of the relevant symmetries of the elliptic force 

feiup (t) = sn [Ktj 7r; m) cn (Ktj n;m), (3) 

where K = K(m) is the complete elliptic integral of the first kind while sn(-;m) and cn(-;m) are Jacobian elliptic 
functions of parameter m 6 [0,1] [16]. Now, m is the single parameter controlling the breakage of the relevant 
symmetries. Unlike functions (1) and (2), function (3) exhibits the advantageous property that its waveform changes 
while its amplitude and image remain constant, f e iu P (t) £ [— 1/2, 1/2] , Vt, as the shape parameter m varies from to 
1 (see Fig. 1 in Ref. [14]). It has been demonstrated [14] that its maximal (with respect to the integration range) 

transmitted impulse over a half-period, I [f] = J T ^ 2 f W dt , and its DSB are monotonously decreasing and increasing 

functions of m, respectively. This result has led to the conclusion that optimal enhancement of DRT is achieved when 
maximal effective (i.e., critical) symmetry breaking occurs, which is in turn a consequence of two reshaping- induced 
competing effects: the increase of the DSB and the decrease of the maximal transmitted impulse over a half-period, 
thus implying the existence of a particular force waveform which optimally enhances DRT. However, reaching the 
same conclusion for the biharmonic force (1) is a harder task due to the aforementioned problematic situation. 




FIG. 1: (Color online), (a) Temporal biharmonic function [Eq. (1)] for the optimal value ifieff ~ tt/2, and (b) the corresponding 
normalized function [Eq. (5)], versus time for different values of r\. 




FIG. 2: (Color online), (a) Temporal biharmonic function [Eq. (1)] for the least favourable value <p e ff = 0, and (b) the 
corresponding normalized function [Eqs. (7) and (8)], versus time t for different values of r\. 

While it has been shown for Hamiltonian systems subjected to symmetric spatial potentials and force (1) [14] that 
when (1 — 77) jr\ ^ 1/2 and <p e ff = {7r/2,37r/2} the DSB is maximal, it is shown in this Letter that for such optimal 
values of tp e ff the maximal transmitted impulse over a half-period exhibits a single maximum at ( 1 — 77) / 77 = 1/2 
(i.e. 77 = 2/3) while the amplitude of the force (1) is held constant, concluding again that optimal enhancement of 
DRT is achieved when critical symmetry breaking occurs. Consider, for example, the case ip e ff — tt/2. Clearly, one 
needs an affine transformation to renormalize the function 77 cos (t) — (1 — rj) sin (2t) in order to change its image from 
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[-M{rj) ,M(r})], where 



3ri + N{r)) 



M(v) = 32 (1 -V) ^ 32 + 30??2 + 2?? [N{Tl) 32] ' 

JVfa) = v/32 - 6477 + 3377 2 , (4) 

to [—1/2, 1/2] for < 77 ^ 1 (in accordance with both the amplitude and the image of the elliptic force (3)). One 
therefore obtains the normalized function 

f* f,\ _ f<Pejf=n/2(t) , . 

W=»/2W- 2M(r?) ' lj 

Figure 1(b) shows a plot of /* =n / 2 (t). One readily finds that its corresponding maximal transmitted impulse over 
a half-period, 

nf; eff =*/2] (v) = j^y, (6) 

exhibits a single maximum at 77 = 2/3, confirming thus the above scenario of optimal enhancement of DRT by critical 
symmetry breaking (see Fig. 4 below for a plot of /[/* ff / 2 ] ~~ ^ versus 77). 

Since the breakage of the shift symmetry of any biharmonic function [Eqs. (1) and (2)] occurs for all values of the 
effective phase, one could expect on the basis of the DSB mechanism that the above scenario should hold for any 
other value of the effective phase, including the least favourable ones, i.e., those values not breaking the corresponding 
relevant reversal symmetry [14]. It is shown here that this is indeed the case by considering the illustrative value 
(p e ff — 0. Using now a different affine transformation, one renormalizes the function r/ cos (t) + (1 — r])cos(2t) to 

change its image from [- ^V-s^ , l] for < r, 4/5 and [1 - 2ry, 1] for 4/5 < r\ < 1 to [-1/2, 1/2] for < J? < 1. 
Noting again that this renormalization process is unique, one straightforwardly obtains the normalized function 




_ j Va „=Q{t)-\/2-R{n)/2 



1 - R( V ) 



(7) 



1 - 277, 4/5 < 77 < 1 J ' 

A plot of f£ =o(t) is shown in Fig. 2(b). One finds that the corresponding maximal transmitted impulse over a 
half-period, 

TC _ ,(,H(I)I±||> „, 

also exhibits a single maximum at r] = 2/3 with I[f^ ff= o\ (v = 2/3) = 7r/6, confirming again the above scenario of 
optimal enhancement of DRT by critical symmetry breaking (see Fig. 4 below for a plot of I[f£ cff=0 ] versus 77). 

If the above scenario is universal (in a sense to be specified below), it must also be found for the spatial force (2) 
since the forces (1) and (2) exhibit the same symmetry properties, and hence the DSB is the same for both types of 
forces. For the sake of comparison with the latter temporal case, consider again the value ip e ff = 0. In this case, 

one has to renormalize the function 77 cos (x) + 2 (1 — 77) cos (2a;) to change its image from 



33^-64i)+32 9 _ 
16-6?7 ' Z 



for 



< 77 8/9 and [2 - 3?y, 2 - 77] for 8/9 77 1 to [-1/2, 1/2] for < 77 < 1. The resulting normalized function is 



now 



(x) = g ^o(.x)-(2-77)/2-P(7 ? )/2 

f 33r, 2 -64r,+32 O < 77 < 8 /9 ] 

2 - 377, 8/9 < 77 1 J 

Figure 3(b) shows a plot of <?* e (x). Accordingly, the corresponding maximal transmitted impulse over a half-period 



eff 

is now 



I K eff =0)(v)={^) 2 _ p( y (12) 
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FIG. 3: (Color online), (a) Spatial biharmonic function [Eq. (2)] for the least favourable value <p e ff = 0, and (b) the 
corresponding normalized function [Eqs. (10) and (11)], versus space x for different values of rj. 




FIG. 4: (Color online). Impulse functions I[f^, Bff =o] (v) [blue (black) line, Eq. (9)], e ,, =0 ] {rj) [red (gray) line, Eq. (12)], 
3.37(/[/* /2 ] - 1) [dashed line, cf. Eq. (6)], and 3.52r7 2 (1 - rj) (dashed-dotted line) versus rj. The inset show the optimal 
waveform cos (z) — 0.5 sin (2z) ,z = {t,x}, for both temporal and spatial biharmonic forces [Eqs. (1) and (2)]. 



which exhibits a single maximum at r\ = 4/5 with J[g* ] (77 = 4/5) = 7r/6, i.e., the same value as 
J[/* e//=0 ] {rj — 2/3) (see Fig. 4 for a plot of /[ffj .. = o] versus 77). Observe that the generality of the present DRT 
scenario is additionally supported by the remarkable conclusion that the two biharmonic forces (1) and (2) have the 
same optimal waveform at their respective optimal values of rj and ip e ff (see Fig. 4, inset). In this sense, such an 
optimal waveform is universal, i.e., once one has identified the relevant spatio-temporal symmetries to be broken in 
a given equation of motion containing the force (1) or the force (2), this optimal waveform maximally enhances the 
ratchet effect in that any other waveform (i.e., any choice of the parameters 77 and f e ff other than the optimal) yields 
a lower ratchet effect while the remaining equation parameters are held constant. 

The present theory explains on the basis of a simple criticality scenario all previously published results for a great 
diversity of systems (see, in particular, the references cited in Ref. [14]). Additionally, this theory is coherent with, and 
explains in a general setting, the general validity of the scaling law {V) ~ ?7 2 (1 — 77) deduced in Refs. [14] and [17] for 
a temporal biharmonic force with small amplitudes. Indeed, one can now expect (V) ~ S {rj) for arbitrary amplitudes, 
where S {rj) exhibits features similar to those of the functions I{f£ eff= o] {rj) and /[<?* ] (77) for temporal and spatial 
biharmonic forces, respectively, i.e., S{rj = 0, 1) = while S {rj) exhibits a single maximum at rj = 2/3 and 77 = 4/5 for 
temporal and spatial biharmonic forces, respectively. It also explains: (i) the effectiveness of the traditionally used 
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ratchet potential V(x) — Vq [sin (2irx/L) + 0.25 sin (4irx/L)] [4-6]; (ii) the experimentally obtained optimal parameters 
of the driving potential, biharmonic in both space and time, of a quantum ratchet [18]; and (hi) the directed ratchet 
transport strength of matter- wave solitons formed in a Bose-Einstein condensate [19]. Experimental confirmation of 
the present findings can be readily obtained, for example in the context of cold atoms in optical lattices [20] . A theory 
providing a universal quantitative criterion to optimize the ratchet effect induced by symmetry breaking of temporal 
and spatial biharmonic forces was demonstrated, paving the way for any future optimal application of the ratchet 
effect. 
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